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{R = La, Sm, Nd, Pr, Ce) 

Gianluca Giovannetti^'^, Sanjeev Kumar^'^, Jeroen van den Brink^'^ 

^Institute Lorentz for Theoretical Physics, Leiden University, 

P.O. Box 9506, 2300 RA Leiden, The Netherlands 

'^Faculty of Science and Technology and MESA+ Research Institute, 

University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands 

^Institute for Molecules and Materials, Radboud Universiteit Nijmegen, 

P.O. Box 9010, 6500 GL Nijmegen, The Netherlands 

(Dated: November 3, 2008) 

We investigate the magnetic phase diagram of the newly discovered iron-based high temperature 
oxypnictide superconductors of the type ROi-^FccFeAs, with rare earths R=La, Sm, Nd, Pr and 
Ce by means of ah initio SGGA and SGGA+U density functional computations. We find undoped 
LaOFeAs to be a Mott insulator when incorporating electronic correlations via SGGA+U for any 
physically relevant value of U . The doped compounds are according to SGGA conductors with 
a transition from an antiferromagnetic to a non-magnetic state at a hole doping of concentration 
Xc=0.075 for R=Nd, Pr and at electron doping Xc=0.25 for Ce and 0.6 for Sm. Superconductivity 
in these rare-earth oxypnictides thus appears in the vicinity of a magnetic quantum critical point 
where electronic correlations are expected to play an important role because of the vicinity of a 
Mott insulating state at zero doping. 
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Introduction Very recently, more than two decades af- 
ter the discovery of high temperature superconductivity 
in the cuprates [1], a new class of iron-based supercon- 
ductors was discovered in the family of rare-earth oxyp- 
nictides [2j |3]. The highest superconducting tempera- 
tures reported so far are for doped oxypnictides of the 
type ROi_^F^FeAs with R=Sm, Nd, Pr and Ce with 
the transition temperatures of Tc =43, 52, 52 and 41 K, 
respectively 01 lU [6l |7] . The Sm compound synthesized 
under high pressure, was recently reported [8 to have a 
Tc of 55 K. The lanthanum compound [2 has a signifi- 
cantly lower Tc of 26 K. 

The electronic structure of these type of oxypnictides 
was first studied with density functional bandstructure 
methods by Lebegue [9], who investigated the related 
phosphide LaOFeP, a material with the same crystal- 
lographic structure as the iron arsenides. Even if this 
compound is build from layers containing formally Fe^+ 
ions, Lebegue showed that this oxypnictide is character- 
ized by large intralayer Fe-P covalency. It was subse- 
quently established that also the electronic bands of the 
arsenide materials are very similar to the ones in the 
phosphide [IQl [TT]. The observed large iron ?>d band- 
width (~4 eV) strongly reduces the tendency for Fe^+ 
to develop a magnetic moment. In the case of LaOFeP 
this leads to the stabilization of a nonmagnetic metal- 
lic groundstate [9 . This result contrasts with recent 
first principles calculations on the arsenide compound 
LaOFeAs. This material was found to be an antiferro- 
magnetic conductor. It was put forward that in the ar- 
senide a spin density wave groundstate forms due to the 
presence of Fermi surface nesting [T2j |13l |14| [15] . Here 
we show that electronic correlations, incorporated within 
SGGA+U drastically change this picture: for physically 
relevant values of U LaOFeAs is a magnetic Mott insu- 
lator. 
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FIG. 1: Schematic view of the crystal structure of ROFeAs 
consisting of layers of corner connected FeAs4 tetrahedra 
(red/yellow) and RO4 tetrahedra (blue/gray). 



The fact that the closely related iron-phosphide and 
arsenide have different magnetic groundstates, suggests 
that at zero temperature these materials could be close 
to an antiferromagnetic - paramagnetic phase transition 
of the iron 3d moments -a quantum critical point. Our 
SGGA ah initio density functional calculations confirm 
that indeed in the rare-earth oxypnictides with a high 
superconducting transition temperature, ROi-a^Fa^FeAs 
with R= Pr, Nd, Sm and Ce, such a quantum critical 
point is present. In the doped compounds supercon- 
ductivity appears in the vicinity of this quantum critical 
point. 



(A) c (A) 



ZR 



ZAs 



SmOFeAs 3.940 8.496 0.1299 0.6402 

NdOFeAs 3.965 8.575 0.1427 0.6437 

PrOFeAs 3.985 8.595 0.1444 0.6417 

CeOFeAs 3.996 8.648 0.1462 0.6426 

LaOFeAs 4.038 8.753 0.1415 0.6512 



TABLE I: Lattice constants and unit cell parameters of RO- 
FeAs [71 125] , with the experimental zr and zas for La 2 and 
values computed by relaxation otherwise. 



I. CRYSTAL STRUCTURE 

The ROFeAs compounds have a ZrCuSiAs type struc- 
ture, characterized by Fe-As layers, where Fe is in the 
center of corner sharing As tetrahedra, see Fig. [5] We 
compute the relaxed atomic positions, the electronic and 
the magnetic structure of these systems within the frame- 
work of the density functional theory [161 HZ] using the 
Vienna ab-initio simulation package (VASP) [T8|. The 
Kohn-Sham equations in the self-consistent calculations 
were solved within a spin generalized gradient approxi- 
mation (SGGA) for the electronic exchange-correlation 
potential (PW91 functional [H]). For the undoped ma- 
terial electronic correlations in the Fe 36Z-shell were ac- 
counted for within SGGA+U [201 HI!- The electronic 
structure is computed using the projector augmented 
wave method (PAW [2?^ '^S^), and the valence pseudo- 
wave-functions are expanded in a plane wave basis set 
with a cutoff energy of 500 eV. All the integrations 
in the Brillouin zone are performed with a tetrahedron 
scheme [24] using a sampling grid of 12x12x6 k-points 
and convergency was checked by performing a number of 
calculations also with a denser grid of k-points. 

The calculations for ROFeAs (R=La, Sm, Nd, Pr, and 
Ce) were performed using the experimental lattice pa- 
rameters [25 . In cases for which experimental values for 
the atomic positions within the unit cell are not avail- 
able, we computed those by minimizing the total en- 
ergy of the paramagnetic state until the relative changes 
in energy were less than 10~^. The combined results 
are summarized in Table L We also computed the re- 
laxed structure in the antiferromagnetic state and within 
the numerical accuracy we find no differences in atomic 
positions. The electronic bandstructures of these sys- 
tems that we calculated are very similar to the one of 
LaOFeAs [lOl HmHllSl [14]. Besides the magnetism 
caused by the iron 3(i moments that is present in all the 
FeAs systems that we have studied, there is another con- 
tribution to the magnetism from the open 4/-shells of 
the rare-earth ions, except for lanthanum. A magnetic 
interplay between those subsystems can certainly occur, 
although probably only at very low temperatures as the 
rare-earth magnetism is generally much weaker than the 
3(i one. Here we wish to focus on the magnetism of the 
Fe-As subsystem and therefore treat the rare-earth 4/ 
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FIG. 2: Gaps A and A^d of LaOFeAs in the magnetic "stripe" 
phase from SGGA+U as a function of U . A is the gap between 
top of the valence and bottom of the conduction band. For 
large values of [/, arsenic bands appear in the gap and closes 
it. Add is the splitting between the occupied and empty Fe 
3d-bands. Also the local iron moment is shown. 



electrons in the computations as core electrons, so that 
they fulfill their role as electron donors but preempting 
their contribution to magnetism. 



II. ELECTRONIC AND MAGNETIC 
STRUCTURE OF LAOFEAS 

For the undoped ROFeAs compounds we find with 
SGGA that R=Nd and Pr are nonmagnetic, while Sm, Ce 
and La are antiferromagnetic, with iron moments of 1.08, 
0.66 and 1.62 /i^, respectively. We also established the 
validity of the premise of earlier computations that the 
related iron phosphide LaOFeP is nonmagnetic [9 . How- 
ever, it is well known that in materials containing transi- 
tion metal ions with open 3(i-shells electronic correlations 
play an important role. We therefore also compute the 
electronic structure for LaOFeAs with SGGA+U. In re- 
cently studied iron compounds the physical relevant val- 
ues of U are between 4 and 6 eV; examples are Fe304 (5 
eV [28 ), FesOBOs (5.5 eV [29]), FeO (4.3 eV, [30;) and 
Fe2Si64 (4.6 and 4.9 eV, [30]). 

The "striped" magnetic state (with ordering vector 
(7r,0) on the squared lattice of iron atoms within the 
FeAs plane) we find in all cases to be the groundstate, 
in agreement with experimental observations [Ml [26l |27] . 
In this magnetic phase already for very small values of U 
the system develops a Mott gap, see Fig. [2] We find for 
instance for U=l eV that A= 33 meV. In the range of 
U between 4 and 6 eV, the Mott gap is around 400-550 
meV. 

A remarkable feature is that the gap closes for large 
values of [/, around 8 eV. This contrast with our finding 
that with increasing U the splitting of the iron d-bands 
increases, as is clear from the value of the splitting be- 
tween the occupied and empty Fe 3(i-bands A^^^^ shown 
in Fig. |2j We obtain /S^dd by computing the dominant 
character of the bands. From this one can estimate /S^^d-, 
but as bands by definition have a mixed character it is 
not rigorously defined. This analysis clearly shows that 
the reason for the closing of the gap A for large U is that 
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FIG. 3: Bandstructure of undoped LaOFeAs in the magnetic 
"stripe" phase, computed with SGGA (left) and SGGA+U, 
using U = 4eV. Bands with predominatly iron 3d character 
are highhghted in black. 



arsenic bands are pushed out of the conduction and va- 
lence bands when U increases. It is evident in the band- 
structure and densities of states that these As bands, 
hybridized with the Fe 3d states, start filling up the d-d 
gap. 

The bandstructure for /7 = and 4 eV is plotted in 
Fig. Is] The bands and the densities of states (DOS) in 
Fig. ^[reveal that the Hubbard U splits the 3d states into 
three sets of bands: a set of bands (I) below 4eV, a set 
of bands, (II) at the top of the valence band and (III) at 
the bottom of the conduction band. We note that there 
is an anomalously large splitting between sets I and III of 
around 6 eV, which is considerably larger than the Hub- 
bard /7 of 4 eV. We will see that this splitting is due to the 
fact that the local Coulomb interactions have two effects: 
{i) they drive a Mott metal-insulator transition and (n) 
for small U an increase of the local moment and there- 
fore significantly increase of local spin splittings. The set 
I corresponds to quasi-localized 3d orbitals, containing 
electrons with parallel spins, giving a large contribution 
to the local Fe moment. The set II corresponds to about 
one electron per Fe with its spin reversed with respect 
to the electrons in the set of bands I. Because of the rel- 
atively large increase of the local moment due to U the 
bands I and II become significantly spin split. The un- 
occupied set of iron bands III is related to empty iron 3d 
states, with a spin direction that is reversed with respect 
to the local moment. These observations are particularly 
evident from the Fe-DOS plotted in Fig. |4] 

For 17=4 eV around the X* point a band is pushed out 
from the valence band that has predominantly arsenic 
character, mixed to a certain extend with iron d. The 
band has a large dispersion and its contribution to the 
DOS is small, see Fig. |4J We observe that for larger 
values of U such a band is pushed out even more until 
for /7 = 8 eV it finally crosses the whole range from 
valence to conduction bands, spanning the energy scale 



> 

^ 

c 

LU 

-2 



-4 



-6 



x2 



r 


} 



x2 



FIG. 4: Density of states (DOS) and projected DOS of un- 
doped LaOFeAs in the magnetic "stripe" phase, computed 
with SGGA (left) and SGGA+U, using U = 4eV. For iron 
the projection along the local Fe magnetic moment is indi- 
cated in black, the one antiparallel to it in red. 



of Add- The detailed origin of this band and its relevance 
for the doped systems is an interesting subject for future 
investigations. 



III. ELECTRONIC AND MAGNETIC 
STRUCTURE OF DOPED IRON-PNICTIDES 

For doped, metallic systems correlation effects can- 
not be included in a meanfield like approach such as 
SGGA+U. For the doped iron-pnictides we therefore only 
use SGGA computations, which are expected to cap- 
ture the effects of Fermi surface instabilities. We find 
that doping changes the ordered moment markedly, see 
Fig.[5| We consider both hole and electron doping, which 
are treated in the computations within the virtual crys- 
tal approximation. We find that upon hole doping the 
antiferromagnetic state appears in the Nd and Pr com- 
pound at a critical concentration Xc=0.075. For both 
the Ce and Sm compound electron doping severely re- 
duces the ordered moment until it vanishes at Xc=0.2b 
and 0.6, respectively, see Fig. [5] The lanthanum system 
also exhibits a reduction of the ordered moment upon 
electron doping, although the antiferromagnetic state is 
stable within the doping range that we consider here. In 
all cases the ferromagnetic state is higher in energy. 

In Fig. |6] the energy difference between the antifer- 
romagnetic and nonmagnetic groundstate is shown to- 
gether with the experimental doping dependence of the 
superconducting Tc of CeOi-a^Fa^FeAs. Clearly the anti- 
ferromagnetic instability extends into the superconduct- 
ing region, with the magnetic quantum critical point cen- 
tered in the superconducting dome. In the other com- 
pounds the critical point is close to the superconducting 
region of the phase diagram. 

It has been pointed out by several authors that the 
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FIG. 5: Ordered iron moment in hole and electron doped 
ROFeAs, computed from first principles. The antiferromag- 
netic (AFM) and paramagnetic (PARA) regions are shown 
for R=Sm, Nd, Pr and Ce. QCP denotes a quantum critical 
point. 
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FIG. 6: Energy difference per unit cell between nonmagnetic 
and antiferromagnetic groundstate of doped CeOFeAs, refer- 
ing to the left vertical scale and the supercoducting transition 
temperatures 7 , refering to the right scale. The magnetic 
quantum critical point at Xc=0.25 is centered around the su- 
perconducting dome. 



antiferromagnetic or spin density wave instability in the 
R=La compound is related to the presence of Fermi sur- 
face nesting between the Fermi surface regions around 
the r and M points in the Brioullin zone [121 1131 IT4[ ITS] . 
The same is found here for R= Sm, Nd, Pr and Ce. The 
doping dependence of this instability is related to the 
sensitivity of the nesting condition to the details of crys- 
tallographic and electronic structure, which we capture 
within a single-particle bandstructure approach. 

It should be noted that electron-electron interactions 
beyond the ones that are effectively taken into account 
in our SGGA computations, are expected to renormalize 
the magnitude of the ordered moments and the related 
critical doping concentrations. Dynamical Mean Field 
Theory results on doped LaOFeAs for moderate values 
of U indicate that electronic correlations do affect the 
hole pockets around the F and electron pockets around 
the M point [11 . 

Another source of renormalization are the transversal 
quantum spin fluctuations which one expects to be par- 
ticularly relevant in these quasi two-dimensional mag- 
netic systems. The treatment of these effects is of course 
beyond the scope of present day density functional the- 
ory and needs to be considered within the framework 
of model Hamiltonians. One might also expect a small 
change in Xc when the impurity potentials of the dopants 
are treated in a scheme beyond the virtual crystal approx- 
imation that we have presented here. 

Nevertheless, our calculations make clearly evident 
that in the high temperature superconducting iron ar- 
senides a magnetic and Mott metal-insulator quantum 
phase transition is present and that superconductivity 
appears in the vicinity of those magnetic quantum critical 
points. The possibility to have a quantum critical state 
of matter close to a magnetic critical point has recently 
attracted much attention because the response of such a 
state is expected to follow universal patterns defined by 
the quantum mechanical nature of the fluctuations ^31^ . 



This universality manifests itself through power-law be- 
haviours of response functions, giving rise to scaling laws 
that are in principle directly accessible experimentally. 



IV. CONCLUSIONS 

We find that undoped LaOFeAs is a magnetic Mott 
insulator with a magnetic "stripe" groundstate for any 
physically relevant value of the Hubbard U within 
SGGA+U, suggesting that the doped materials are close 
to a zero temperature metal-insulator transition. This 
conclusion is supported by a recent analyses of electronic 
transport data by Si and Abrahams |32j. The Mott gap 
depends on the precise value of U and is estimated to be 
between 400 and 550 meV. In addition our SGGA calcu- 
lations show that doping the rare-earth oxypnictides RO- 
FeAs with R= Sm, Nd, Pr and Ce induces a phase transi- 
tion from an antiferromagnetic to non-magnetic ground- 
state. 

It is a priori not clear whether the magnetic instabil- 
ity competes with the superconducting one and thus sup- 
presses superconductivity or whether it is the root cause 
of the superconducting state. The fact that it was con- 
cluded on the basis of calculations of the electron-phonon 
coupling strength that a phonon mediated mechanism 
for superconductivity is unlikely to be operative in these 
compounds [15, 33 (although contested by others [34] ) 
suggests that in these materials the fluctuations associ- 
ated with the magnetic quantum-critical point are essen- 
tial for superconductivity. Such a relationship between 
magnetism, quantum criticality and superconductivity 
has already been observed in a number of heavy fermion 
materials, for instance CeRhlns [35 , CeCu2Si2 j36l 157] 
and UPd2Al3 [38 and uranium intermetallic UCoGe [39] . 
This suggests a close relation of these materials to the su- 
perconducting rare-earth oxypnictides. 
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